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antisite  -  mercury  vacancy  pair  in  HgCdTe.  We  find  a  binding  energy 
of  more  than  1  eV. 
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1.  NATIVE  POINT  DEFECTS  IN  HgCdTe  AND  RELATED  IR  MATERIALS 

We  have  calculated  the  binding  energy  for  the  tellurium  antisite  mercury  vacancy 
pair,  corresponding  to  the  reaction 

+  TeHg  (VHgTeHg)po,y .  (1) 

where  V}ig  and  Teng  are  the  isolated  mercury  vacancy  and  tellurium  antisite  point 
defects  and  (VHgTeHg)po,v  is  the  bound  mercury  vacancy  tellurium  antisite  pair  in 
which  the  point  defects  reside  on  adjacent  cation  sites.  We  expected  this  pair  to 
be  important  for  two  reasons.  First,  both  the  mercury  vacancy  and  the  tellurium 
antisite  are  important  defects  in  Hgo.7sCdo.22Te  and  as  isolated  species  are  present  in 
significant  quantities.  Second,  the  mercury  vacancy  is  an  acceptor  and  the  tellurium 
antisite  is  a  donor,  and  therefore  a  Coulomb  attraction  will  contribute  to  their  binding. 
Finally,  the  tellurium  antisite  produces  a  compressive  strain  and  the  mercury  vacancy 
produces  a  tensile  strain  (although  somewhat  smaller  in  magnitude  than  that  of 
the  antisite),  and  therefore  there  will  be  a  mechanical  attraction  between  the  point 
defects. 

The  binding  energy  of  the  pair  was  calculated  using  the  self-consistent  full  po¬ 
tential  linearized  muffin-tin  orbital  (FP-LMTO)  method  and  supercells  discussed  in 
previous  quarterly  reports.  We  have  completed  the  calculation  using  the  32-atom 
supercell  for  both  the  isolated  d^ects  and  the  defect  pair  and  with  a  minimal  basis 
set,  and  we  obtained  a  binding  energy  of  1.1  eV.  We  are  currently  completing  the 
calculations  with  a  larger  basis  set  size  and  with  a  larger  supercell.  These  calculations 
are  extremely  computationally  intensive,  and  strategies  to  expedite  the  calculations 
are  being  developed. 

We  have  generalized  our  previously  developed  quasichemical  method  to  include 
defect  complexes.  We  are  currently  completing  a  preliminary  analysis  of  the  density 
of  this  defect  pair  we  expect  in  Hgo.7BCdo.22Te.  We  are  also  completing  the  analysis  of 
this  defect  in  light  of  the  fact  that  the  tellurium  antisite  is  likely  to  be  a  slow  diffuser, 
and  therefore  may  not  reach  equilibrium  either  as  an  isolated  defect  or  as  a  part  of  a 
defect  complex.  Defects  of  this  analysis  will  be  presented  in  the  next  quarterly  report. 
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2.  CALCULATION  OF  IONIZATION  ENERGIES 


An  important  part  of  our  calculation  of  defect  densities  in  semiconductors  (HgCdTe, 
ZnSe,  HgZnTe,  and  CdTe)  is  the  determination  of  the  ionization  energies.  In  HgCdTe, 
we  used  a  rather  crude  analysis  to  determine  whether  the  native  point  defects  were 
donors  or  acceptors  because  our  calculations  of  the  densities  of  the  defects  were  rela¬ 
tively  insensitive  to  the  exact  determination  of  the  ionization  energies.  Unfortunately, 
this  crude  analysis  does  not  yield  the  exact  position  of  the  levels  in  the  gap  for  either 
the  first  or  second  ionization  levels,  nor  does  it  indicate  whether  a  defect  is  a  negative- 
U  center.  For  the  wider-gap  semiconductors  such  as  CdTe  and  ZnSe,  the  locations 
of  the  levels  in  the  gap  are  much  more  important  and  will  significantly  impact  the 
densities  of  the  defects.  Thus,  we  are  working  towards  a  method  to  better  locate  the 
defect  levels  in  the  gap.  Much  of  the  work  this  quarter  has  been  focused  on  this  task, 
and  is  summarized  below. 

We  assume  that  the  total  energy  of  a  solid  is  a  smooth  fvmction  of  the  system 
charge,  z,  that  is 


E{Z)  =  Eq  -f-  E!iz  •+  E^z^ .  (2) 

The  energy  per  electron  to  change  the  system  charge  from  zi  to  Z2  is  to  lowest  order 
given  by 

z.)  =  =  £,  +  E,{z,  -  »,)•  (3) 

Z2  —  Zi 

If  zi  and  zj  are  integers  differing  by  unity,  Ei^^  is  just  the  ionization  energy.  Note 
that 

(It) 

Now  for  an  infinite  system,  the  change  in  energy  when  an  electron  is  added  or  removed 
from  the  system  is  the  Fermi  energy  or  electron  chemical  potential  tp  so  that 

(5) 

We  are  interested  in  the  one-electron  ionization  levels  of  defects,  that  is,  the 
energy  of  the  levels  with  respect  to  the  band  edges.  We  will  always  take  the  energies 
with  respect  to  the  valence  band  edge  due  to  the  local  density  approximation  (LDA) 
band-gap  errors  (consequences  of  which  on  the  defect  levels  are  to  be  discussed  further 
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below).  For  the  one-electron  levels  we  take  the  difference  between  the  Fermi  level  for 
the  half-filled  state  and  the  energy  of  top  of  the  vedence  band,  fvi  that  is 


Ei(l)  =  €f(z  =  —0.5)  —  Cv 

(6) 

Ed(2)  =  cf(z  =  —1.5)  —  Cv 

for  the  first  and  second  ionization  levels  of  a  donor  and 

(7) 

Bail)  =  €e  -  Cf(*  =  +0.5) 

(8) 

Ea{2)  =  Cc  “  €f(«  =  +1.5) 

(9) 

for  the  first  and  second  ionization  levels  of  an  acceptor.  Here  ep  is  the  zero-temperature 
Fermi  level  of  the  system  with  a  nondegenerate  density  of  a  single  defect  species,  is 
taken  from  the  bulk  calculation,  and  Cc  is  determined  from  -I-  Because 

the  zero  of  energy  is  arbitrary  in  the  LDA  calculations,  we  have  assured  a  proper 
alignment  of  the  bands  (and  consequently  the  Fermi  energies)  by  aligning  bulk-like 
characteristic  features  in  the  density  of  states  of  the  bulk  and  defect  calculations. 

We  have  corrected  for  the  finite  width  of  the  gap  state  in  the  following  manner. 
We  know  that  in  the  limit  of  very  large  systems  the  gap  states  will  be  very  narrow, 
approaching  a  delta  function,  and  will  accept  two  electrons  per  defect  for  s-like  (or 
Al-symmetry)  states  and  six  electrons  for  p-like  (or  T2-synunetry)  states.  Because 
we  are  doing  supercell  calculations,  our  localized  levels  in  the  gap  are  quite  broad 
because  of  the  overlap  of  the  wave  functions  of  point  defects  in  adjacent  supercells. 
We  correct  for  this  finite  width  by  associating  the  ionization  energy  level  with  the 
center  of  gravity  of  this  localized  level  (as  determined  from  the  density  of  states). 

2.1  NATIVE  POINT  DEFECTS  IN  CdTe 

We  applied  the  analysis  discussed  above  to  determine  the  energy  levels  of  native 
point  defect  in  CdTe.  Table  1  presents  our  preliminary  results  for  CdTe. 

A  preliminary  analysis  of  the  density  of  the  native  point  defects,  including  both 
neutral  and  ionized  states,  has  been  completed.  We  find  that  in  equilibrium,  the 
cadmium  vacancy  and  the  tellurium  antisite  dominate  at  low  cadmium  pressures, 
while  the  cadmium  vacancy  and  the  cadmium  interstitial  dominate  at  the  higher 
cadmium  pressures.  A  more  complete  analysis  will  be  presented  in  the  next  quarterly 
report. 


Table  1.  Ionization  energies  for  the  native  point  defects  in  CdTe.  Results  are 
from  LMTO  calculations  using  the  atomic-spheres  approximation,  and  do  not  include 
Jahn-Teller  distortions. 


Defect 

Acceptor  levels 

Donor  levels 

single 

J2B23I 

Vcd 

0. 

0.10 

Vt. 

0.70 

0.85 

Cdxe 

1.30 

1.40 

Tecd 

isa 

0.33 

Cd/cd 

1.00 

Cd/.p, 

iSi 

0.56 

Te/cd 

in 

1.20 

0.55 

0.84 

2.2  NATIVE  POINT  DEFECTS  IN  ZnSe 

We  are  still  completing  the  calculation  of  the  ionization  energies  of  the  native 
point  defect  in  ZnSe.  Although  we  have  not  yet  completed  the  ionization  energy 
calculation,  we  have  demonstrated  that  equilibrium  native  defects  can  compensate 
the  p-type  dopants  in  ZnSe. 

Figure  1  shows  the  neutral  native  point  defects  in  ZnSe  at  600  as  a  fimction 
of  zinc  partial  pressure.  The  zinc  vacancy  is  the  dominant  defect  followed  by  the 
seleniiuu  vacancy  at  high  zinc  pressures  and  the  selenium  antisite  at  low  zinc  pressures. 
From  our  calciilations  to  date,  the  zinc  vacancy  appears  to  be  an  acceptor,  while  both 
the  selenium  antisite  and  vacancy  appear  to  be  donors. 

To  demonstrate  the  p-type  compensation  mechanism,  we  consider  two  cases.  First 
we  assume  that  the  zinc  vacancy  has  a  single  acceptor  level  at  midgap,  while  the  sele¬ 
nium  vacancy  and  antisites  are  assumed  to  be  double  donors  with  their  first  ionization 
level  at  the  conduction  band  edge,  and  the  second  level  at  0.1  eV.  Using  these  ioniza¬ 
tion  levels  and  the  defect  formation  energies  from  our  completed  calculations  for  the 
neutral  defects  (used  to  generate  Figure  1),  we  have  completed  a  defect  analysis  at 
600  with  just  these  three  defects.  We  find  that  the  material  can  be  easily  doped 
n-type  throughout  the  phase  stability  region,  but  a  high  level  of  p*type  doping  cannot 
be  achieved.  This  situation  is  illustrated  in  Figure  2,  where  we  have  completed  our 
calculations  of  the  defect  concentrations  for  an  impurity  concentration  of  10^°  cm— 3. 


4 


zinc  partial  pressure  (atm) 


(a)  with  10  ■•9  cm*3  acceptor  impurities 


Figure  2  Total  native  point  defect  concentrations  (neutral  and  ionized)  in  ZnSe,  which  is  equilibrated  at 
600°C  assuming  Vse  and  Sezn  are  shallow  donors  and  Vzn  is  a  deep  acceptor.  Electronic  equilibration 
at  room  temperature 


Defect  concentration  (cm-3) 


Next,  we  have  considered  the  case  in  which  the  zinc  vacancy  has  associated  with 
it  two  acceptor  levels  at  0  and  0.1  eV,  while  the  selenium  antisite  and  vacancy  are 
assumed  to  have  singly  ionizable  midgap  levels.  Results  for  p-  and  n-type  impurity 
concentrations  of  10^^  cm~^  are  illustrated  in  Figure  3.  As  one  can  see  in  this  case, 
the  material  is  easily  doped  p*type  but  cannot  be  doped  n-type. 

Obviously  the  presence  of  self- compensation  in  ZnSe  depends  on  the  true  position 
of  the  ionization  levels  in  the  gap.  Our  preliminary  results  indicate  that  neither  of 
the  extreme  limits  discussed  above  will  be  valid.  Our  present  results  indicate  that  the 
acceptor  level  associated  with  the  zinc  vacancy  will  be  relatively  shaUow,  and  because 
of  the  dominance  of  the  neutral  zinc  vacancy  (Figure  1),  we  expect  that  this  defect 
will  be  the  major  player  in  the  compensation  game.  Thus,  at  this  point  we  expect 
that  the  native  point  defects  are  not  responsible  for  the  self-compensation  mechanism 
in  ZnSe. 

When  our  native  point  defect  analysis  is  complete,  we  will  begin  to  examine  issues 
related  to  the  dark  line  defects  in  ZnSe,  their  relationship  to  the  stacking  faults,  and 
means  for  their  elimination  and/or  passivation. 

3.  HgZnTe  AS  AN  LWIR  DETECTOR  MATERIAL 

The  transport-related  properties  such  as  electron  mobility.  Hall  coefficient,  and 
Fermi  levels  are  calctdated  with  accurate  analytical  band  structures,  Feimi- Dirac  (FD) 
statistics  and  experimental  energy  gap.  Thus,  calcdated  values  differ  substantially 
from  the  ones  obtained  with  effective  mass  or  parabolic  band  structme  approximation 
of  Hgo.rsZno.irTe  alloy.  We  find  that  approximating  the  Hall  factor  by  unity  over 
a  wide  range  of  carrier  concentrations  is  accurate  only  for  temperatures  less  than 
100  K.  An  error  of  about  50%  is  expected  at  higher  temperatures  and  low  carrier 
concentration  (10^*  cm"®),  and  about  20%  is  expected  at  low  T  and  high  carrier 
concentration  (10^*  cm"®).  The  impurity,  phonon  and  alloy  disorder-limited  mobility 
calculations  with  full  solution  to  Boltzmann  transport  equation  with  FD  statistics 
predict  a  hump  near  40  K.  The  mobility  in  HgZnTe  is  about  10-20%  higher  than  that 
in  HgCdTe. 

3.1  INTRODUaiON 

The  experimental  results  on  electron  transport  properties  of  semiconductors  is 
often  compared  to  theory  that  rests  on  three  approximations,  namely,  parabolic  band 
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zinc  partial  pressure  (atm) 

(a)  with  10'*^  cm'3  acceptor  impurities 


Figure  3  Total  native  point  defect  concentrations  (neutral  and  ionized)  ZnSe  which  is  equilibrated  at  600° 
assuming  Vse  and  Sezn  are  deep  donors  and  Vzn  is  a  shallow  acceptor.  Electronic  equilibration  at  room 
temperature 


Defect  concentration  (cm-3) 


structures  for  those  states  occupied  in  the  measurement,  Maxwell-Boltzmann  (MB) 
statistics,  and  collision  time  approximation  to  the  full  Boltzmann  gain-loss  equation. 
These  approximations  are  made  to  all  scattering  mechemisms  whether  they  are  elastic 
or  inelastic.  It  is  well  known  that,  even  in  large-gap  materials,  the  constant  effec¬ 
tive  mass  approximation  is  valid  only  very  near  (»  EgjlQ)  to  the  band  edge  (Kane, 
1957;  Krishnamurthy  et  al.,  1987).  This  approximation  has  been  recognized  to  be 
particularly  poor  for  narrow-gap  materials,  and  nonparabolic  corrections  calculated 
in  the  k.p  formalism  are  often  used  (Schmidt  1970,  Meyer  and  Bartoli,  1982;  Bartoli 
et  al.,  1982).  This  correction  is  substantial  but  still  differs  considerably  from  our 
more  accurately  calculated  band  structures.  In  addition,  our  fit  of  the  conduction 
band  to  an  analytical  function  makes  many  results  transparent  and  simplifies  the 
calculations.  As  the  Fermi  energy  can  easily  move  into  the  conduction  band  of  lightly 
doped  small-gap  materials,  the  form  of  the  Boltzmann  equation  with  FD  (instead  of 
the  usual  MB)  statistics  must  be  used  to  obtain  accurate  transport  coefficients. 

Here  we  report  restdts  from  our  study  of  absorption  coefficient,  Fermi  energy. 
Hall  coefficient,  and  electron  mobility  calculated  with  FD  statistics  and  accurate 
pseudopotential  band  structures  fine  tuned  with  tight-binding  (TB)  corrections. 

3.2  BAND  STRUCTURE 

Quantitatively  accurate  band  structures  of  group  IV  elements  (Krishnamurthy  et 
al.,  1986),  III-V  compounds  (Chen  and  Sher,  1981),  and  II- VI  compoimds  (Chen  and 
Sher,  1981;  Berding  et  al.,  1987)  can  be  obtained  using  a  minimum  set  of  sp^  orbitals 
in  semiempirical  calculations.  First,  empirical  pseudopotential  form  factors  are  used 
to  calculate  a  TB  Hamiltonian  H  in  the  minimum  set.  H  is  then  transformed  into  a 
zeroth-order  Hq  in  an  orthonormal  basis.  Then,  a  perturbative  Hamiltonian  having  a 
first-neighbor  TB  form  is  added  to  Hq  to  fine  tune  the  band  structure.  Because  long- 
range  interactions  are  included  in  this  Hamiltonian,  the  measured  band  curvatures 
are  correctly  reproduced.  This  procedure  is  followed  for  both  HgTe  and  CdTe,  and 
then  the  alloy  band  structures  are  calculated  in  the  coherent  potential  approximation. 

We  focus  Hgo.22Cdo.78Te  alloy  with  100-meV  band  gap  for  the  studies  reported 
hare.  We  find  that  the  calculated  conduction  band  is  replicated  very  well  by  a  hyper¬ 
bola. 


Et  =  {ik’’ +  <?)''’- c 


(10) 


11 


where  7  and  c  are  adjusted  to  fit  the  calculated  band  structure  in  the  energy  range  of 
interest.  When  7  and  c  are  treated  as  constants  related  to  the  band  gap  Eg  and  the 
effective  mass,  this  expression  will  reduce  to  the  same  nonparabolic  correction  form 
obtained  in  the  k-p  method.  Our  values  of  7  and  c  are  47.78  and  0.058,  respectively, 
whereas  the  corresponding  k*p  values  are  36.0  and  0.05.  The  differences  are  found 
to  be  large  enough  to  cause  a  noticeable  change  in  the  band  structure.  The  band 
structure  calculated  by  diagonalizing  the  Hamiltonian  is  shown  in  Figure  1  (thick 
line).  We  find  that  our  hyperbola  fits  the  band  structure  calculated  by  diagonalizing 
the  HaTniltonia.il  quite  well  up  to  an  energy  of  0.5  eV  firom  the  conduction  band  edge. 
Without  loss  of  accuracy,  in  the  studies  considered  here,  Eq.  (10)  is  used  as  the 
energy  dispersion  relation  in  transport  expressions  that  follow. 

For  energies  E  —  Ec  greater  than  50  meV  where  the  shape  of  the  conduction  band 
is  nearly  linear  in  k,  the  group  velocity  is  a  constant  independent  of  the  k.  Then 
the  density  of  states  (DOS)  increases  proportional  to  E  rather  than  E^^^  as  in  the 
case  of  parabolic  bands.  Clearly,  these  features  modify  the  transport  properties  of 
electrons  occupying  these  states.  As  we  will  show  in  the  following  section,  at  the 
carrier  concentration  and  temperatures  often  found  in  device  structures,  the  Fermi 
level  falls  into  the  region  where  these  features  contribute. 

3.3  FERMI  LEVEL 

The  calculation  of  Fermi  level  ejras  a  function  temperature  T  and  doping  concen¬ 
tration  njD  is  required  for  all  transport  calculations.  A  knowledge  of  temperature- 
dependent  gap  Ef(T)  is  essential  to  obtain  cf.  Ideally,  temperature  dependence 
should  be  included  in  the  Hamiltonian  from  which  the  variation  of  E,  with  T  could 
be  obtained.  We  have  developed  a  general  method  to  incorporate  phonon  effects  on 
the  band  gap.  However,  we  have  used  the  experimental  values  (Sher  et  al.,  1986)  of 
the  gap.  After  studying  the  effects  of  various  approximations  and  obtaining  trends, 
the  calculations  will  be  repeated  with  temperature-dependent  Hamiltonian. 

Then,  the  tpis  calculated  from  the  condition  (Sze,  1981)  that  at  a  given  T  the 
number  electrons  in  the  conduction  band  is  the  sum  of  electrons  excited  from  the 
valence  band  and  donor  levels.  In  this  study  where  the  modifications  caused  by  the 
band  structures  are  being  emphasized,  the  donor  states  are  assumed  to  be  located 
at  the  bottom  of  the  conduction  band.  The  valence  and  conduction  band  DOS  are 
calculated  from  our  band  structure.  The  valence  band  DOS  yields  a  hole  effective 
mass  of  0.65.  The  eF(measured  from  valence  band  edge)  as  a  function  of  T  and  ud 


calciilated  from  hyperbolic  band  structures  are  given  in  Figure  4.  Hyperbolic  band 
makes  substantial  change  from  parabolic  band  generated  values. 

3.4  HALL  COEFFICIENT 

The  carrier  density,  n  in  n-type  material  is  normally  deduced  from  measurements 
of  Hall  coefficient  given  by  re/en,  by  assuming  the  Hall  factor  is  unity.  If 
one  uses  parabolic  approximation  and  MB  statistics,  is  approximately  unity.  We 
set  out  to  examining  the  effect  of  removing  these  approximations  using  the  correct 
band  structures.  It  requires  generalizing  the  Boltzmann  transport  equation  (BTE)  to 
include  FD  distribution  function.  We  have 

^  £  l«,(k,  k')/(k')(i  -  /(k))  -  tt(k',  k)/(k)(i  -  /(k-))]  (11) 

The  first  term  of  the  right-hand  side  is  the  gain  term,  and  the  second  term  is  the  loss 
term.  In  equilibrium,  the  left-hand  side  of  Eq.  (11)  is  identically  zero  and  f  becomes 
the  equilibrium  FD  distribution  function  fo,  given  by 


/o(£») (12) 

where  ^  is  (k^T)**^.  In  the  presence  of  electric  and  magnetic  field, 


dm  _  dm 

di  dt 


+  ^f{k)-^(E  +  vxB) 


(13) 


In  steady  state,  the  d((k)/dt  in  Eq.  (13)  vanishes.  In  the  small  field  regime,  we  can 
linearize  f  and  write  it  as  a  sum  of  fb  and  a  perturbation  fi(k).  Disregarding  the 
derivative  of  fi(k)  and  after  some  algebraic  manipulation,  Eq.  (13)  reduces  to. 


V/o(k)  •  j(E  +  V  X  B)  =  E  [W(k,  k')/.(k')  -  W(k',  k)/i(k)l  (14) 

where  the  renormalized  W  and  the  usual  transition  probability  per  unit  time  w  are 
related  by 


W(k,k')  =  «.(k,k')^-|M  (15) 

Note  that  for  elastic  scattering  W  and  w  are  equal.  However,  for  inelastic  cases 
the  effect  depends  on  whether  energies  at  k  and  k’  are  larger  or  smaller  than  ep.  If 
both  initial  and  final  energies  are  larger  (or  smaller)  than  ep,  only  a  small  correction 
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Figure  4  Fermi  energy  as  a  function  of  temperature:  HgZnTe 


to  w  is  expected.  However,  if  the  initial  state  is  above  e/rand  the  final  state  is  below 
6jr,  that  scattering  is  suppressed. 

In  the  collision  time  approximation,  the  gain  term  in  Eq.  (15)  is  neglected  and 
the  effective  collision  time  is 


(^)-'  =  EW(kMc)  (16) 

k' 

Using  this  collision  time  approximation,  it  is  straightforward  to  obtain  the  expression 
for  r*: 


fe  =  ZkgT 


ri:k/ci:i.*M(rf)V.(i-/o) 


(17) 


where  VkJS*  is  7k  k.  The  calculated  r*  with  hyperbolic  band  structures  is  shown 
for  various  T  and  doping  densities  in  Figure  5.  The  r#  is  approximately  1  for  lower 
temperatures.  However,  at  higher  temperatures,  it  shows  considerable  structure  with 
a  maximum  value  of  about  1.6.  The  ratio  of  re  calculated  with  hyperbolic  bands  to 
that  obtained  with  parabolic  bands  is  shown  in  Figure  6.  Our  calculations  show  that 
parabolic  values  can  be  in  error  by  as  much  as  50%.  The  experimentalists  normally 
use  a  value  of  1  for  ig  in  deducing  the  intrinsic  density  &om  Hall  measurements, 
suggesting  that  reported  values  may  be  smaller  by  as  much  as  50%  at  intermediate 
temperatures. 

3.5  DRIFT  MOBILITY 


The  formalism  developed  above  to  find  the  solution  to  BTE  with  FD  statistics  can 
be  used  to  calculate  the  mobility  fi.  We  obtain  the  following  expression  in  collision 
time  approximation: 


ZrkeT  [ 


El.*W/o(l-/o) 


(18) 


The  calculated  mobility  by  including  the  scattering  due  to  ionized  impurities  and 
polar  optic  phonons  as  a  function  of  temperature  and  doping  density,  is  shown  in 
Figure  7.  The  change  in  the  Debye  screening  length,  complicated  variation  of  Fermi 
level,  and  phonon  scattering  give  rise  to  a  crossover  in  mobility  near  40  K,  as  seen 
in  Figure  7.  However,  it  is  well  known  that  the  collision  time  approximation  always 
overestimates  the  velocity  transition  rate  that  determines  the  mobility.  For  mean* 
ingful  comparison  with  experiment,  we  generalized  the  above  procedure  to  get  full 
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solution  to  BTE  (Krishnamurthy  and  Sher,  1993).  Thus,  the  calculated  mobility  with 
our  hyperbolic  band  structure  and  usual  k‘p  band  structure  is  shown  along  with  that 
obtained  in  collision  time  approximation  and  compared  with  values  in  HgCdTe  in 
Figure  7. 

One  important  feature  of  Figure  7  is  that  the  hump  in  mobility  near  40  K  could  be 
explained  with  competing  impurity  and  phonon  scattering  rates.  A  hump  in  mobility 
is  HgCdTe  is  seen  experimentally.  Also  note  that  the  mobility  in  HgZnTe  is  about 
10-20%  higher  than  in  HgCdTe. 

3.6  CONCLUSIONS 

We  have  studied  the  effect  of  various  approximations  on  electron  transport  co¬ 
efficients  and  on  ways  to  extract  physical  parameters  &om  experiments.  We  point 
out  how  the  values  interpreted  from  experiments  depend  crucially  on  vwious  approx¬ 
imations  such  as  effective  mass,  MB  statistics,  and  collision  time.  The  main  results 
are 

•  Approximating  the  Hall  factor  by  unity  over  a  wide  range  of  carrier  concentrations 
is  accurate  only  for  low  temperatures,  typically  less  than  100  K.  An  error  of  about 
50%  is  expected  at  high  temperattires  and  low  carrier  concentration  (10^®  cm"®),  and 
about  20%  is  expected  at  low  T  and  high  carrier  concentration  (10^*  cm"®). 

•  The  mobility  calc\ilation8  with  fvill  solution  to  BTE  with  FD  statistics  predict 
a  hump  near  40  K.  The  mobility  in  HgZnTe  is  about  10-20%  higher  than  that  in 
HgCdTe. 


4.  WORK  PLANNED 

In  the  next  quarter  we  will  continue  the  calculation  of  the  ionization  energies  of 
the  defects  in  CdTe,  ZnSe,  and  HgCdTe  and  the  incorporation  of  these  energies  into 
the  quasichemical  formalism  to  predict  densities  of  the  ionized  defect  concentrations. 
We  will  continue  the  refinement  of  the  calculation  of  the  tellurium  antisite  mercury 
vacancy  defect  pair.  We  will  present  several  papers  at  the  MCT  Workshop  in  October, 
and  will  be  preparing  for  these  presentations. 
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